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From December 2011 to May 2014, about 5 fb“^ of data were taken with the BESIII detector at center- 
of-mass energies between 3.810 GeV and 4.600 GeV to study the charmoniumlike states and higher excited 
charmonium states. The time-integrated luminosity of the collected data sample is measured to a precision of 
1% by analyzing events produced by the large-angle Bhabha scattering process. 

PACS numbers: 13.66.Jn 


I. INTRODUCTION 


As a r-charm factory, the BESIII experiment has collected 
the world’s largest sample of e+e” collision data at center-of- 
mass (CM) energies between 3.810 GeV and 4.600 GeV. In 


this energy region, the charmoniumlike states and higher ex¬ 
cited charmonium states are produced copiously, which makes 
comprehensive studies possible. 

The charmoniumlike states discovered in recent years have 
drawn great attention of both theorists and experimentalists 



3 


for their exotic properties, as reviewed e.g. in Ref. m. Being 
well above the open charm threshold, the strong coupling of 
these states to hidden charm processes makes their interpre¬ 
tation as conventional charmonium states very difficult. On 
the other hand, the theory of the strong interaction. Quantum 
Chromodynamics (QCD), does not prohibit the existence of 
exotic states beyond the quark model, e.g. molecular states, 
tetraquark states, hybrid states, etc. Either the verification or 
the exclusion of the existence of such states will help to evalu¬ 
ate the quark model and better understand QCD. Even though 
some states have been identified as higher excited charmo¬ 
nium states, such as the 'tp{4:0A0), '0(4160), and 0(4415), their 
large widths and the interference with each other make their 
precise study complicated. In addition, the relationship be¬ 
tween the charmoniumlike states and higher excited charmo¬ 
nium states is still not clear. The precise knowledge of the 
time-integrated luminosity is essential for quantitative analy¬ 
sis of these states. 

In this paper, we present a measurement of the integrated 
luminosity based on the analysis of the Bhabha scattering pro¬ 
cess e+e” —>■ ( 7 )e“''e“. A similar method has been used in 
the luminosity measurement of 0(3770) data at BESIII 12 • 
The process has a simple and clean signature and a large pro¬ 
duction cross section, which allows for a small systematic and 
a negligible statistical uncertainty. A cross check of the result 
is performed by analyzing the di-gamma process e+e” —^ 77 . 


II. THE DETECTOR 


BESIII is a general purpose detector which covers 93% of 
the solid angle and operates at the e+e” collider BEPCII. A 
detailed description of the facilities is given in Ref. 0. The 
detector consists of four main components: (a) A small-cell, 
helium-based main drift chamber (MDC) with 43 layers pro¬ 
vides an average single-hit resolution of 135 /rm, and a mo¬ 
mentum resolution of 0.5% for charged tracks at 1 GeV/c 
in a 1 T magnetic field, (b) An electro-magnetic calorime¬ 
ter (EMC), consisting of 6240 CsI(Tl) crystals in a cylindrical 
structure (barrel and two endcaps). The energy resolution for 
1.0 GeV photons is 2.5% (5%) in the barrel (endcaps), while 
the position resolution is 6 mm (9 mm) in the barrel (end- 
caps). (c) A time-of-fight system (TOE), constructed of 5 cm 
thick plastic scintillators, arranged in 88 detectors of 2.4 m 
length in two layers in the barrel and 96 fan-shaped detec¬ 
tors in the endcaps. The barrel (endcap) time resolution of 
80 ps (110 ps) provides laKj'n: separation for momenta up to 
about 1.0 GeV/c. (d) A muon counter (MUC), consisting of 
nine layers of resistive plate chambers in the barrel and eight 
layers for each endcap. It is incorporated in the iron return 
yoke of the superconducting magnet. Its position resolution is 
about 2 cm. A GEANT4 lIHISl based detector simulation pack¬ 
age has been developed to model the detector response. Due 
to the crossing angle of the beams at the interaction point, the 
e+e” CM system is slightly boosted with respect to the labo¬ 
ratory frame. 


III. DATA SAMPLE AND MONTE CARLO SIMULATION 

Twenty-one data samples have been taken at CM energies 
between 3.810 GeV and 4.600 GeV. Six of the data sets exceed 
the others in accumulated statistics by an order of magnitude. 
These samples were taken on the peaks of charmoniumlike 
states, like the Y(4260), Y(4360), and Y(4630), or higher ex¬ 
cited charmonium states, like 0(4040), and 0(4415), in order 
to study these resonances and their decays in great detail. The 
data samples taken at the other CM energies serve as scan 
points to study the behavior of the cross section around these 
resonances. All individual data samples are listed in Table|^ 

At each energy point, one million Bhabha events were gen¬ 
erated using the BABAYAGA3.5 i) generator with the options 
presented in Table [H] Eor the BABAYAGA3.5 generator, the 
uncertainty in calculating the cross section is 0.5% , which 
meets the demand of the total uncertainty of luminosity mea¬ 
surement. The kinematic distributions of the final state par¬ 
ticles from the BABAYAGA3.5 generator are consistent with 
those from data. In the simulation, the scattering angles of the 
final state particles were limited to a range from 20 degrees to 
160 degrees, which slightly exceeds the sensitive volume of 
the detector, in order to save on computing resources. An en¬ 
ergy threshold of 0.04 GeV was applied on the final state par¬ 
ticles. The acolinearity of the events has not been constrained. 
Einally, the generation was taking into account the running of 
the electromagnetic coupling constant and final state radia¬ 
tion (ESR). 

To study the background and optimize the event selection 
criteria, an inclusive Monte Carlo (MC) sample corresponding 
to a luminosity of 500 pb“^ at the CM energy of 4.260 GeV 
was generated, in which the QED processes, the continuum 
production of hadrons, and the initial state radiation (ISR) 
to J/0 and 0(3686) resonance process were included. The 
BABAYAGA3.5 generator was used to simulate the QED pro¬ 
cesses, signal and background. Other processes, such as the 
decays of the J/0, were generated with specialized models 
that have been packaged and customized for the BESIII Of¬ 
fline Software System (BOSS) (see 0 for an overview). 


IV. EVENT SELECTION AND RESULT 

Signal candidates are required to have exactly two oppo¬ 
sitely charged tracks. The tracks must originate from a cylin¬ 
drical volume, centered around the interaction point, which is 
defined by a radius of 1 cm perpendicular to the beam axis 
and a length of ±10 cm along the beam axis. In addition, the 
charged tracks are required to be within | cos 9\ < 0.8, where 
9 is the polar angle, measured by the MDC. Without applying 
further particle identification, the tracks are assigned as elec¬ 
tron and positron depending on their charge. The deposited 
energies of electron and positron in EMC must be larger than 
X 1.55 (GeV) to remove the di-muon background, where 
^/s is the CM energy in GeV; the momenta of electron and 
positron are required to be larger than x 2 (GeV/c), to 
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TABLE I. Center-of-mass energy, luminosity obtained from the nominal measurement (L), cross check results (Lck), and relative differences 
between the two results. The uncertainties are statistical only. Superscripts indicate separate samples acquired at the same CM energy. 


CM energy (GeV) L (pb 

^ck (pb Relative difference (%) 

3.810 

50.54±0.03 

50.11±0.08 

-0.85±0.17 

3.900 

52.61±0.03 

52.57±0.08 

-0.08±0.17 

4.009 

481.96±0.01 

480.54±0.23 

-0.30±0.05 

4.090 

52.63±0.03 

52.37±0.08 

-0.49±0.17 

4.190 

43.09±0.03 

43.08±0.08 

-0.03±0.20 

4.210 

54.55±0.03 

54.27±0.09 

-0.62±0.18 

4.220 

54.13±0.03 

54.22±0.09 

+0.17±0.18 

4.230^ 

44.40±0.03 

44.64±0.08 

+0.54±0.20 

4.230^ 

1047.34±0.14 

1041.56±0.37 

-0.56±0.04 

4.245 

55.59±0.04 

55.52±0.09 

-0.13±0.18 

4.260^ 

523.74±0.10 

524.57±0.26 

+0.16±0.06 

4.260^ 

301.93±0.08 

301.11±0.20 

-0.28±0.08 

4.310 

44.90±0.03 

45.29±0.08 

+0.87±0.19 

4.360 

539.84±0.10 

541.38±0.28 

+0.29±0.06 

4.390 

55.18±0.04 

55.27±0.09 

+0.16±0.18 

4.420^ 

44.67±0.03 

44.77±0.08 

+0.22±0.20 

4.420^ 

1028.89±0.13 

1029.63±0.37 

+0.07±0.04 

4.470 

109.94±0.04 

109.51±0.13 

-0.39±0.13 

4.530 

109.98±0.04 

109.47±0.13 

-0.46±0.13 

4.575 

47.67±0.03 

47.57±0.08 

-0.21±0.18 

4.600 

566.93±0.11 

563.45±0.28 

-0.62±0.06 


TABLE 11. Options for the BABAYAGA3.5 generator used to gener¬ 
ate the simulated MC data samples. 


Parameters 

Value 

Ebeam 

2.130 GeV or others 

MinThetaAngle 

20 ° 

MaxThetaAngle 

160° 

MinimumEnergy 

0.04 GeV 

MaximumAcollinearity 

180° 

RunningAlpha 

1 

ESR switch 

1 


suppress background events from lighter vector resonances 
produced in the ISR process, such as J/ij), V'(3686) and other 
resonances, decaying into e+e“ pairs. For the data sample 
with a CM energy of 3.810 or 3.910 GeV, the effect of the re¬ 
maining '0(3686) events is studied by applying a 20% larger 
momentum requirement, and is found to be negligible. The 
requirements on the deposited energies and momenta are not 
optimized in detail, as the number of the signal events in such 
an analysis is large enough. All the variables mentioned above 
are determined in the initial e+e” CM frame. The ratio of the 
number of remaining background events to the number of sig¬ 
nal events, estimated from the inclusive MC sample, is found 
to be less than 2 x 10“^, which is negligible. Thus all the 
selected events are taken as Bhabha events. 

Figure[2shows the comparisons between data and MC sim¬ 
ulation for the kinematic variables of the leptons by taking 
data at the CM energy of 4.260 GeV as an example. Rea¬ 
sonable agreement is observed in the angular and momentum 


distributions. The striking difference between data and simu¬ 
lation found in the distributions of energies deposited by the 
leptons in the EMC emerges from imperfections in the simu¬ 
lation of the energy response of individual detector channels. 
At the CM energies analyzed in this work, a single shower 
in the calorimeter can be so energetic that the deposited en¬ 
ergy per crystal exceeds the dynamic range of the analog-to- 
digital converter (ADC), causing individual ADC channels to 
saturate. In the analysis presented in this report, conditions 
applied on the energy deposits are not affected. Relevant de¬ 
viations between data and MC are considered as contributions 
to the systematic uncertainties. 


The integrated luminosity is calculated with 

ATohs 

Bhabha 
f^Bhabha ^ ^ 

where A^g^ltha is the number of observed Bhabha events, 
CTBhabha is the cross section of the Bhabha process, and e is 
the efficiency determined by analyzing the signal MC sample. 
The cross sections are calculated with the BABAYAGA3.5 gen¬ 
erator using the parameters listed in Table[n]and decrease with 
increasing energies. The efficiencies are almost independent 
of the CM energy, as intended by the choice of relative con¬ 
ditions on lepton momenta and deposited energies. The lumi¬ 
nosity results calculated with Equation are listed in Table 
The statistical accuracy of the resulting integrated luminosity 
is better than 0.1% at all energy points. 
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FIG. 1. Comparison between data and MC simulation at the CM energy of 4.260 GeV. The top row is for positron and the bottom row for 
electron. From left to right, the plots show the distribution of deposited energy in EMC, the distribution of the cosine of the polar angle 
measured hy the MDC, and the distribution of the track momentum from the MDC. Black dots with error bars illustrate data and the red ones 
are MC simulation. Note that the y-axis is in logarithmic scale and the MC is normalized to data by the number of events for each sub-plot. 
When drawing the distribution of one variable, the requirements on the other variables are applied. 


V. SYSTEMATIC UNCERTAINTY 


The following sources of systematic uncertainties are con¬ 
sidered; the uncertainty of the tracking efficiency, the uncer¬ 
tainty related to the requirements on the kinematic variables, 
the statistical uncertainty of the MC sample, the uncertainty of 
the beam energy measurement, the uncertainty of the trigger 
efficiency, and the systematic uncertainty of the event genera¬ 
tor. 

To estimate the systematic uncertainty related to the track¬ 
ing efficiency, the Bhabha event sample is selected using in¬ 
formation from the EMC only, without using the tracking in¬ 
formation in the MDC. The selection criteria are: at least two 
clusters in the EMC for each candidate, and the two most ener¬ 
getic clusters are assumed to originate from the pair; the 
deposited energies of the two clusters are required to be larger 
than X 1.8 (GeV). At CM energies above 4.420 GeV, 

fs 

the requirement is changed to x 1.55 (GeV). This ad¬ 
justment allows to avoid additional systematic uncertainties 
which would be introduced by the deviation of data and sim¬ 
ulation in the deposited energy in the EMC, as discussed in 


Sec.IV. The polar angle of each cluster is required to be within 
|cos0™®| < 0.8, where polar angle measured 

by the EMC; to remove the background from the di-photon 
process, Ai/) is required to be in the range of [—40°, —5°] or 
[5°, 40°], where A(/) = \(j)i — (f >2 \ — 180° and (j)i ^2 are the az¬ 
imuthal angles of the clusters in the EMC boosted to the CM 
frame. The efficiency that the selected Bhabha events pass 
through the track requirements applied in the nominal anal¬ 
ysis is calculated for both data and MC sample, and the dif¬ 
ference between them is taken as the systematic uncertainty 
connected to the tracking efficiency. 

The systematic uncertainty in the requirement on the po¬ 
lar angle is estimated by changing the requirement from 
|cos0| < 0.8 to |cos0| < 0.7. The difference between the re¬ 
sulting luminosity and nominal one is taken as the associated 
systematic uncertainty. The systematic uncertainty caused by 
the requirement on the energy deposited in the EMC is esti- 

/s 

mated by changing the requirement from x 1.55 (GeV) to 

fs 

X 1.71 (GeV). The systematic uncertainty caused by the 
requirement on the momentum is estimated by changing the 
requirement from x 2 (GeV/c) to x 2.06 (GeV/c). 
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The ranges are picked as these cause the largest deviations 
from the nominal luminosity result near the requirements ap¬ 
plied. 

The statistical uncertainty of the efficiency determined from 
MC simulations is 0.25%. The CM energy is determined us¬ 
ing —>■ events. The invariant mass of the di¬ 

muon system is calculated taking into account ISR and FSR 
effects 0 . The difference between the CM energy listed in 
Table and the one measured with di-muon process is about 
2 MeV, and the corresponding systematic uncertainty is esti¬ 
mated by changing the CM energy by 2 MeV in the MC simu¬ 
lation. The trigger efficiency for the Bhabha process is 100% 
with an uncertainty of less than 0 . 1 % 0. The theoretical un¬ 
certainty of the cross section calculated by the BABAYAGA3.5 
generator is given as 0.5% ( 6 ]. 


The same systematic uncertainty estimation method is ap¬ 
plied to all the sub-samples. The largest relative uncertainty 
among them is taken as the associated uncertainty for all the 
sub-samples. The systematic uncertainties considered in this 
work are summarized in Table III By assuming the sources of 
the systematic uncertainties to be uncorrelated, the total un¬ 
certainty is calculated as 0.97% by adding the contributions in 
quadrature. 


TABLE III. Summary of the systematic uncertainties. 


Source 

Relative uncertainty (%) 

Tracking efficiency 

0.39 

Energy requirement 

0.09 

Momentum requirement 

0.43 

Polar angle requirement 

0.38 

MC statistics 

0.25 

Beam energy 

0.42 

Trigger efficiency 

0.10 

Generator 

0.50 

Total 

0.97 


VI. CROSS CHECK 


To verify the result, a cross check with di-gamma events is 
performed. The event selection criteria are the same as those 
used in estimating the systematic uncertainty caused by the 
tracking efficiency, except for the requirement on A(j). In or¬ 
der to reduce the Bhabha background, the Acj) is required to be 
in the range of [— 0 . 8 °, 0 . 8 °], since photons are not deflected 


in the magnetic field. 

The luminosity results of this cross check (Lck) shown 
in Table together with the relative differences to the nom¬ 
inal ones. Both results are well consistent for all individual 
measurements, indicating the robustness of the result. 

VII. SUMMARY 


The integrated luminosity of the data samples taken at BE- 
SIII for studying the charmoniumlike states and higher ex¬ 
cited charmonium states is measured to a precision of 1 % with 
Bhabha events. The total uncertainty is dominated by the sys¬ 
tematic uncertainty. A cross check with di-gamma events is 
performed and the results are consistent with each other. The 
result presented here is essential for future measurements of 
cross sections with these data, and it has already been used in 
the discovery of charged charmoniumlike states IIoHSl. 
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